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Mass loss and expansion of ultra compact dwarf galaxies through 
gas expulsion and stellar evolution for top-heavy stellar initial mass 
functions 
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ABSTRACT 

The dynamical U-band mass-to-light ratios of ultra compact dwarf galaxies (UCDs) are higher 
than predicted by simple stellar population models with the canonical stellar initial mass func- 
tion (IMF). One way to explain this finding is a top-heavy IMF, so that the unseen mass is 
provided by additional remnants of high-mass stars. A possible explanation for why the IMF 
in UCDs could be top-heavy while this is not the case in less massive stellar systems is that 
encounters between proto-stars and stars become probable in forming massive systems. How- 
ever, the required number of additional stellar remnants proves to be rather high, which raises 
the question of how their progenitors would affect the early evolution of a UCD. We have 
therefore calculated the first 200 Myr of the evolution of the UCDs, using the particle-mesh 
code Superbox. It is assumed that the stellar populations of UCDs were created in an initial 
starburst, which implies heavy mass loss during the following » 40 Myr due to primordial 
gas expulsion and supernova explosions. This mass loss is modelled by reducing the mass of 
the particles according to tabulated mass loss histories which account for different IMFs, star 
formation efficiencies (SFEs), heating efficiencies (HEs), initial masses and initial extensions 
of the computed UCDs. For each combination of SFE and HE we find objects that roughly 
resemble UCDs at the end of the simulation. For low SFEs, the IMF would have to be steeper 
than in the case of very high SFEs for the models not to expand too much. However, the main 
conclusion is that the existence of UCDs does not contradict the notion that their stellar pop- 
ulations formed rapidly and with a top-heavy IMF. We find tentative evidence that the UCDs 
may have had densities as high as 10 8 M Q pc -3 at birth. This will have to be confirmed by 
follow-up modelling. 
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1 INTRODUCTION 

Ultra compact dwarf galaxies (UCDs) are stellar systems with total 



radii of < 50 pc (Hilker et al|l999 


; Drinkwater et alJl2000Ll2003l: 


PhilliDDS et alj2001; Hasesan et al. 


2005|). Thev can be considered 



to be galaxies because of their high median two-body relaxation 
times, t r h, which are at least of the order of a Hubble time, th, 
whi le star clusters, including globular cluste rs (GCs), have t T h < 
th dKroupalll998l : lDabringhausen et alj|2008l) . 

One of the most intriguing prope rties of UCDs are their gen- 
erally high dynam ical M/Lv ratios ( Dabringhause n et al. 2008; 
Mieske et al. 2008). Different explanations have been suggested 
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for this finding, such as the presence of non-ba ryo nic cold dark 
matte r (CDM) in them (e.g. lHasegan et alJliooll and lGoerdt et al.l 
2008) o r the disturbance of UCD s by the tidal field of a massive 
galaxy l lFellhauer & K roupa 2006). However, if dwarf spheroidal 
galaxies (dSphs) are indeed DM dominatecQ and if UCDs are lo- 
cated at the centre of the same type of haloes as dSphs, the DM- 



1 There is an ongoing debate on the origin of the dSphs around the Milky 
Way. Their disk-like distribution has a natural explanation if the dSphs are 
ancient tida l dwarf galaxies instead of DM-dominated primordial galaxies 
iMetz et al. 2009 and references therein). The high M / L ratios derived for 
them would in this scenario either be the consequence o f the assumption 
of virial equilibrium not holding for them I Kroupa 1997) or would imply 
that Newtonian gravity cannot be applied in the limit of very weak fields. 
A tidal origin of dSphs may suggest the same for dwarf elliptical galaxies, 
since [kormendy 1 1985) argues that these two populations may actually be 
the same type of galaxies. 
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density in UCDs would be two orders of magnitude too low to 
explai n their elevated M/L y ratios, although adiabatic contrac- 
tion telumenthaletalJll986h may alleviate this problem jMurravl 
2009). Tidal distortion can explain the high M/Ly ratios of only 
a few UCDs out of a larger sample, as it requires quite specific 
orbital parameters in order to have an observable effect. On the 
other hand, the massive star cluster W3 in the merger remnant 
galaxy NGC 7252 has a mass and a projected half-light radius typ- 
ical for a UCD, while its age suggests th at it formed during th e 
merger of the progenitors of NGC 7252 jMaraston et alj|2004h . 
Fellha uer & Kroupal ^2005) have shown that star cluster complexes 
as observed in interacting sytems like the Antennae (NGC 4038 
and NGC 4039) are likely to evolve into an object similar to W3 on 
the required time-scale, but stellar systems orig inating from tidal 
intera ctions would essentially be CDM-free dBarnes & Hernquist 
In summary, an unusual stellar initial mass function (IMF) 
appears to be an attractive and physically plausible alternative for 
explaining the M / L ratios of UCDs. 

The IMF is a function defining the mass spectrum of stars 
born in a single star-formation event. If age, metallicity and IMF 
of a stellar population are known, its M/Lv ratio can be calcu- 
lated. For a given metallicity and a high enough age, a high M/Ly 
ratio of a stellar population would either indicate an IMF with 
very many low-mass stars (bottom-heavy IMF) or an IMF with 
very many high-mass stars (top-heavy IMF). In the case of a top- 
heavy IMF, the high M/Ly of the stellar population is the con- 
sequence of a high number of stellar remnants, which contribute 
mass, but almost no V-band luminosity. As an explanation for the 
M/Ly ratios in UCDs, a b ottom-heavy IMF has been discussed 
in lMieske & Kroupal 1 20081) . while a top-heavy IMF has been dis- 
cussed in lMurravl ( l2009^ and Dabringhausen, Rroupa & Baumgardt 
(2009), hereafter DKB. 

Proposing a variability of the IMF might seem daring at first 
sight, because so far surveys of st ars have failed in providing 
supportive evidenc e for this notion I Rroupa 2001; Kroupal [20021 ; 
iKumar et al.l 120081) . This finding implies an invariant, universal 
IMF, which is referred to as the canonical IMF. It can be written 



£ c (m„) = kknn, ai , 
with 

1.3, fei = 1, 



(1) 



Q2 



2.3, k 2 = fci 0.5° 



0.5, 



0.1 «: 
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< 0.5, 

^ 77l ma x , 



where m* is the initial stellar mass, m max is the upper mass limit 
of the IMF, the factors hi ensure that the IMF is continuou s where 
the p ower changes and k is a normalisation constant (Kroupa 
2008). The subscript c identifies the canonical IMF. £, c (ni*) equals 
if m, < 0.1 Mq or m* > m max . For stellar systems as massive 
as the UCDs, m max is equal to the maximum mass for stars, which 
is close to 1 50 Mq JWeidner & Kroupall2004l ; IOev & Clarkel2005l ; 
lFigerll2005h . For any IMF, dN = £(?n*)drn* is the number of 
born stars in the mass interval [m*,m* + dm*]. In the present 
paper, the constant k is chosen such that 



5(m,)ra, dm* = 1 
Using this normalisation, 
V = / £(m*)dm* 



(2) 



(3) 



tipliying equation [3] by the factor M^q/Mq therefore equals the 
initial number of stars in a star cluster with an initial stellar mass of 
M*,o and the mean stellar mass, m, equals equation l(2j divided by 
equation l[3}. 

Note that there are limitations to the determination of the IMF 
from star counts. For instance, massive stars are short-lived, which 
is why this approach can only give the high-mass IMF for recent 
star formation events. Low-mass stars on the other hand can be 
almost as old as the Universe, but they can only be detected very 
locally. 

The existence of a universal law for the stellar mass spectrum 
would indeed be surprising from a theoretical point of view, since 
models for star-formation predict that the stellar mass spectrum de- 
pends on the conditions under which star formation takes place 
(e.g. 



; on the conditions under which star tormation takes place 
Adams & Fatuzzol 1 19961 iMurrav & Linl 1 199(1 iLarsonl 1 19981 



and IClarketal .1120071) . Moreover, a top-heavy IMF is in fact re 
quired in a number of astrophysical models. This includes, besides 
the model proposed in DKB for the UCDs , also models for globular 
clusters (GCs) (D 'Antona & Caloil 120041 ; IPrantzo s & Cha rbonneJ 
20061; iDecressin et alj|2007lfl distant galax i es dBaugh et alj|2005l ; 
Nagashima et al 120051; Ivan DokkurrJ 120081; ICharvl l2008lT arid the 



is formally the number of stars whose total mass is 1 Mq . Mul- 



Galactic centre l Maness et al.ll2007F) 7 The motivations for the top- 
heaviness of the IMF in these models include a higher ambient tem- 
perature at the time when the observed population formed and vio- 
lent star formation in particularly dense gas. These conditions were 
likely to be given in the young UCDs, since the universe was much 
younger when they formed (i.e. the temperature of the cosmic mi- 
crowave background was higher) . Furthermore, the a-enrichment 
found bv lEvstigneeva et ahl ( 120071) in most of the Virgo-UCDs sug- 
gests rapid star-formation. 

However, if stellar remnants are to account for the unseen 
mass in the UCDs, the top-heavyness of the IMF would have to 
be very pronounced. Introducing an IMF that equals the canonical 
IMF below 1 M but has a different slope, a, for m > 1 Mq, DKB 
suggest 1.0 < a < 1.6, depending on the age of the UCDs. These 
high-mass IMF slopes imply that the clear majority of the total ini- 
tial stellar mass was locked up in stars more massive than 8 Mq, in 
contrast to the case with the canonical IMF. These stars have a very 
high luminosity and evolve rapidly, which makes their abundance 
a key issue for the evolution of a stellar system. 

If UCDs indeed are the most massive star clusters, their stellar 
populations would essenti ally have formed in a single burst over a 
time-span of ~ 1 Myr (cf. lElmegreenl2 000 ; Hart mann et alj200lh . 
meaning that their stars evolve almost simultaneously. Considering 
the high energies involved in massive star evolution, this implies 
that UCDs with very top-heavy IMFs (with high-mass IMF-slopes 
1.0 < a < 1.6, see DKB) could have lost 90% of their initial stel- 
lar mass over a time span of ~ 40 Myr (which is the lifetime of the 
least massi ve stars that evolve into SNe, cf. the stellar evolution- 
ary grid by |Schalleretai]|l992h . If there was residual gas (i.e. gas 
that was not used up in star formation) in them, which was swept 
out during this phase of violent star cluster evolution, the mass loss 
would have been even more pronounced. Such an extensive mass 
loss shapes the later appearance of a stellar system and may even 
be critical for its survival, if it happens on a short enough time scale 



2 If UCDs are indeed the most massive GCs , as proposed for in stance in 
iMieske et al] l2002l) , lMieske et alj 12004 and lForbes et"ai] 120081) , it is ev- 
ident that a top-heavy IMF in GCs suggests the same for UCDs. Note how- 
ever that residual gas expulsion fr om mass-segregated c lusters alleviates the 
need of a top-heavy IMF in GCs I Decressin et al. 2008). 
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Table 1. The initial param eters and some derived quantities for the UCD-models. The last two lines show the models for newly formed ONC-type star clusters 
from lKroupa et alj 1200 ll) for comparison. The columns denote the identification number of the model, the initial Plummer-radius -R p i,o, its total initial mass 
Mpi o, its stellar initial mass M» t o, the star formation efficiency (SFE= Mo,*/Af p i q)> the heating efficiency, the initial characteristic crossing-time T cr , the 
initial characteristic three-dimensional velocity dispersion <X3rj o> the initial central mass density and an estimate for the time-scale on which a given proto-star 
encounters a star during the formation of the UCD (see Section |2.1.2> . 



model 


^pl,0 


M p i,o 




M t ,o 




SFE 


HE 


I cr 


°"3D,0 


Ppl,0,c 


ienc 




[pc] 


[M ] 




[M ] 








[Myr] 
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m7_r3_sl_hl 


3.0 


1.0 x 
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1.0 x 
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1.0 


1.0 


0.153 


65.0 


0.088 


0.23 


m7_r5_sl_hl 


5.0 


1.0 x 


10 7 


1.0 x 


10 7 


1.0 


1.0 


0.330 


50.3 


0.019 


1.4 


m8_r5_sl_hl 


5.0 


1.0 x 


10 x 


1.0 x 


10 s 


1.0 


1.0 


0.104 


159.0 


0.191 


0.043 


m7j3-s04_hl 


3.0 


2.5 x 


10 7 


1.0 x 


10 7 


0.4 


1.0 


0.097 


102.8 


0.221 


0.15 


m7.r5.s04.hl 


5.0 


2.5 x 


10 7 


1.0 x 


10 7 


0.4 


1.0 


0.209 


79.6 


0.048 


0.87 


m8.r5_s04_hl 


5.0 


2.5 x 


10 s 


1.0 x 


10 s 


0.4 


1.0 


0.066 


251.6 


0.477 


0.027 


m7_r3_s04_h003 


3.0 


2.5 x 


10 7 


1.0 x 


10 7 


0.4 


0.03 


0.097 


102.8 


0.221 


0.15 


m 7_r5_s04_h003 


5.0 


2.5 x 


10 7 


1.0 x 


10 7 
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0.03 


0.209 


79.6 


0.048 


0.87 
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5.0 


2.5 x 


10 s 


1.0 x 


10 s 


0.4 


0.03 


0.066 


251.6 


0.477 


0.027 


ONC A 


0.345 


1.12 x 10 4 


3.75 x 10 3 


0.33 




0.23 


6.8 


0.065 


9.4 


ONCB 


0.158 


1.25 ) 


< 10 4 


4.17 ; 


< 10 3 


0.33 




0.066 


10.8 


0.759 


0.52 



jBoilv & KrouDal2003l : lFellhauer & Kroupal2005h . However, DKB 
argued from structural parameters that mass loss on a time scale of 
40 Myr for UCDs is probably in the adiabatic regime and therefore 
inflates them, but does not threaten to dissolve them. It is clear that 
a numerical study of this issue, including a more detailed treatment 
of mass loss through stellar evolution and residual gas expulsion, is 
necessary to confirm these arguments. It is provided in this paper. 



2 SETUP 

2.1 Initial conditions 

In the present paper, UCDs are assumed to have formed in the 
monolithic collapse of a fragmenting gas clou d, and thus in contrast 
to the mod el for UCD-formation propo sed in lFellhauer & Kroupal 
d2002l) and Fell hauer & Kroupal J2005t) . i.e. the merger of a star 
cluster complex into a single object (see also Kroupa 1998). This 
is not to say that the merging of star clusters is completely irrele- 
vant for UCD formation. For instance, the densest part at the cen- 
tre a proto-UCD could undergo monolithic collapse, while in the 
outskirts of the proto-UCD a multitude of star clusters is formed, 
which eventually merge. However, the apparent universality of the 
IMF in star clusters below the mass-scale of a UCD suggests that 
a UCD could not have a different IMF, if it is exclusively build up 
from such systems. 

The adopted formation scenario for UCDs thus suggests that 
they are the most massive star clusters, which implies that their 
stellar population form ed rapidly in ~ 1 Myr (cf. Elme greet j 2000l : 
iHartmann et al. 1200 lh. The a-enrichment of the UCDs in Virgo re- 
ported by lEvstigneeva et al] j2007h indeed suggests a short time 
scale for star formation, although for UCDs in other environments, 
this q-enrichment is less pronounced or even absent (cf. fig. 8 in 
iMieske et al.l2007l) . For simplicity, we assume that the stellar popu- 
lations of UCDs have formed instantaneously instead of over a very 
short time-span. This is a conservative assumption for the present 
study, since it focusses on the stability of UCDs. A stellar popu- 
lation that is built up over an extended time-span also releases the 
total energy it produces (through stellar processes) over a longer 
time-span. In consequence, the mass loss from UCDs, which is 



powered by the energy produced by the stellar population, will be 
slower and therefore less threatening for the stability of the UCD. 



2.1.1 Structural parameters 

The UCD-models are set up with their ma ss distributed accord- 
ing to the Plummer-model jPlumme3ll91lh . The Plummer-model 
is the simplest plausible and self-consistent mod el for a star cluster 
dBinnev & Tremainel [l987: Heggie & Hut 2003). The advantage of 
using Plummer-models is that all major quantities are analytically 
accessible. 

We choose nine different combinations of initial stellar mass, 
M*,o, initial Plummer-radius, i? p i,o, star formations efficiencies 
(SFEs) and heating efficiencies (HEs) for the UCD-models. The 
choices of the mentioned parameters are detailed below and the 
considered combinations of them are listed in TableQ] together with 
some major quantities derived from them. 

M tj o is chosen such in the models that a stellar mass of the 
order of 10 6 to 10 7 M remains after the evolution of the mas- 
sive stars has come to an end. This is the range in which the stel- 
lar masses of the observed UCDs lie. The chosen values for 72 p i,o 
are either 3 pc or 5 pc and th us similar to the o bserved radii of 
GCs (eg. lMcLaugMmll2000l or ljordan et aHl2005t) . This leads to 
initial central densities, p p i,o,c> ranging from 1.9 x 1O 4 M pc -3 
to 4.8 x 1O 5 M pc -3 (Table[T|l. These values for p p i,o,c are simi- 
lar to the ones that have been calculated for Galactic open clusters, 
such as the Orion Nebula Clus t er (ON C), whose initial parameters 
are discussed in iKroupa et al] ( 1200 lh . They consider models with 
P P i,o,c = 6.5 x 1O 4 M pc" 3 or p pl , , c = 7.6 x 1O 5 M pc" 3 for 
that star cluster, as can be calculated from the initial masses and 
half-light radii given in their table 1. The models discussed here 
are thus not extreme because of the densities in their central re- 
gions, but because of the extension of this central region. This may 
account for the proposed top-heaviness of the IMF in UCDs, see 
Section |2~l~2l 

Embedded star clusters in the Milky way are thus less 
extended than the models discussed in this paper (also see 
lLada & Ladj2003l their table 1). Note that also GCs were initially 
less extended than the UCD-models discussed in this paper, un- 
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less they lost very little mass since their formation. The reason why 
smaller i? p i,i are not considered here are the extreme initial central 
densities they would imply for the objects (for instance of the or- 
der of 10 7 Mq pc -3 for Rpifi = 1 pc; also see Fig. I14t. Besides, 
the very small crossing times of such objects would make compu- 
tations of the evolution very time-consuming while mass loss from 
them would approach the adiabatic regime, where the behaviour of 
the cluster can be calculated analytically with equation below. 

The actual values of the SFE and the HE are hard to quantify. 
In order to get an idea of how these parameters would influence 
the early evolution of a UCD, vastly different and in some cases 
extreme values for them are considered in this paper. 

The SFE is defined as the fraction of the primordial gas that is 
converted into stars during a star-forming event within the cluster- 
or UCD-forming cloud core region. In the UCD-models, it is taken 
to be 1 or 0.4, the latter value being approxi mately the upper limi t 
of the SFEs reported for open star clusters jLada & Ladall2003t) . 
These high choices for the SFEs in UCDs are motivated by the 
fact that it would be more difficult to expel the primordial gas from 
UCDs than from o pen clusters because of the de ep potential wells 
of UCDs (see also El megreen & Efremovlll997|). It has even been 
suggested (e.g. in Elmegreenl l 19991 orfMurrav 2009) that all avail- 
able gas is turned into stars, if the forming stellar system is dense 
and massive enough. If indeed all star clusters and UCDs form on 
the same time scale, the star formation rate must be higher in UCDs 
than in any of the less massive stellar systems. Taking 1 Myr as the 
characteristic time scale for star formation in these systems, the 
average star formation rate in UCDs would be 10-100 M0yr -1 
(DKB). Assuming UCDs are essentially star clusters and that star 
formation is the more rapid the denser the primordial gas is, this 
could be understood if the primordial gas cloud forming a UCD 
is, compared to open star clusters, compressed to a higher density 
during its collapse. 

The HE is defined as the fraction of the energy released by 
stellar processes that actually drives gas out of a star-forming region 
instead of being radiated away. That is, the HE is the ratio between 
the kinetic energy of the interstellar medium (ISM) expelled from 
the stellar system to the total energy inserted into its ISM. The HEs 
in starbursts have been argued to be near 1 in some studies (e.g. 
IChevalier&Clegdll985h , while others suggest that only a few per- 
cent of the energy inserted into the IS M is turned in kinet ic energy 
of gas leaving the stellar system (e.g . iRecchi et al.|[200ll ; also see 
iMelioli & de Gouveia Dal Pinoll2004r) . For the present paper, HEs 
of 1 and 0.03 are considered. 

A major improvement compared to the rather arbitrary choice 
of SFEs and HEs made here would clearly be to estimate these 
parameters in self-consistent modelling of a collapsing gas cloud 
large enough to form a UCD. This would also clarify how long it 
would actually take the stellar population to form in such a system, 
but is currently not a computable option. 



2.7.2 A possible influence of encounters on the IMF 

If UCDs indeed formed with the initial conditions proposed here, 
the likeliness for close encounters between members of their 
emerging stellar populations (stars and proto-stars) would be what 
sets them apart from ONC-like star clusters. This motivates why 
the IMF in UCDs might be top-heavy, while this is not observed in 
star clusters like the ONC. 

The case of a proto-star encountering a star is of particular 
interest. A proto-star exists over a time o f 10 5 yr until most of 
its mass has accreted onto the central core JWuchterl & Tscharnuted 
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H — ' 

a> 
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central density [ solar masses / pc ] 

Figure 1. The expected time until a proto-star collides with a star, t enc , 
in a forming UCD, assuming that half of its total stellar population has 
already been formed (see Section l2.1.2l for details). The estimated t en c for 
the UCD-models, in which all gas is eventually converted into stars, are 
indicated by the three (blue) points. The grey shaded area is where t cnc is 
below the approximate life-time of a proto-star, which is assumed to be 10 5 
years. The solid lines show the t cnc as a function of p p i,o,c for different 
constant M p i o, starting from 10 3 Mq and increasing by a factor of 10 
downwards. The dotted lines show the t cnc as a function of p p i,o,c f° r 
different constant R p ifi, with Rpifi being 0.1, 0.3, 0.5, 1, 3 and 5 pc from 
top to bottom. 



120031) and is thus short-lived, compared to the characterisic time- 
scale for star-formation in a star cluster (~ 1 Myr). It has however 
a radius, r prot o, of the order of 100 AU for essentially all stellar 
masses, since the dependency o f r pro t on the mass of the proto-star 
is only weak (cf. equation 4 in I Goodwin et al 1 120071) . This makes 
an encounter of a proto-star with a star quite likely, as soon as a 
considerable stellar population is already present. 

To estimate a characteristic time-scale for such an encounter 
for the UCD-models listed in Table Q] consider Plummer-spheres 
with the initial parameters from that table. Their density-profiles 



are given as 



p(R) = 



3M p i,o 



1 + 



R 

Rp\,o 



(4) 



(equation 8.51 in lKroup3l2008h . and thus their central stellar den- 
sities at the time when half of their stellar populations have formed 
can be estimated as 

3M„,o 



Ppl.*,c — 



(5) 



This density implies a volume that contains one star on average, 
V». It can be written as 



V, = 



(6) 



where m is the average stellar mass. The time-dependent volume 
through which a proto-star in the central region has travelled due to 
its motion can be written as 
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Figure 2. As Fig. (T}, but for the UCD-models, where a star formation ef- 
ficiency of 0.4 instead of 1 is assumed. While UCD-models shown here 
are assumed to have the same stellar masses as the ones shown in Fig. {T} 
their total masses are higher, leading to higher velocity dispersions and thus 
shorter t cnc at the same stellar density. The meaning of the solid and the 
dotted lines is the same as in Fig. {T}, but the lowest constant value for 
Mpi o chosen here is 2.5 X 10 3 Mq and increasing by a factor of 10 down- 
wards with every solid line. The open squares sh ow models A and B for the 
initial states of ONC-type star clusters from iKroupa et al l fcOOll) . which 
have, compared to the UCD-models in this figure, a slightly lower SFE of 
0.33. 



V(t) = 7rr proto (J3D,o,c t, 

where (73d,o.c is the central 3D velocity dispersion, which is 



C3D,0,c 



GM p i,o 



(7) 



(8) 



2i? P i,o 

where G is the gravitational constant (cf. equation 8.59 in[Kroupa 
2008). The values calculated from equation {8]l for the UCD- 
models in this paper are an order of magnitude higher than in 
the models for the initial states of ONC-type star clusters from 
IKroupa etailfcOOlh . while their central densities are essentially the 
same (see Table [TJ. The time t = t cnc by which a proto-star is to 
be expected to have encountered a star can be calculated be setting 
V(t) — V, and solving for t. Thus, 



r pr°to Ppl,*,c 03D,O,c 



(9) 



Assuming that a top-heavy IMF results from an canonical IMF by 
the collisions of proto-stars with stars, m = 0.65 Mq (which is the 
average stellar mass for a canonical IMF, see Table[2j, and r pro to = 
100 AU = 4.85 x 10 -4 pc are reasonable choices for emerging 
open star clusters and emerging UCDs alike. This implies that the 
IMF would be canonical until it is altered under the influence of 
encounters between the members of an emerging stellar population 
and would stay canonical in stellar systems where such encounters 
are rare at all times. 

Note that the derivation of equation l[9) implicitly assumes 
that the cross section for an encounter of a proto-star with a star is 



the geometrical cross section, A g 



whereas the actual 



gravity. If both the proto-star and the star have the same mass, m, 
this actual cross section is given as 



A = 7rr proto (l + 9) 



(10) 



at the centre of the emerging UCD, where O is the Safronov num- 
ber, 



9 = 



2G1 



(11) 



cross section for such an encounter is higher due to the influence of 



°"3D,0,c r proto 

( lMurrav&Lmi ri996). However, assuming r pro to = 100 AU and 
m = O.65M0 leads to = 0.26 in the less compact model for the 
initial states of ONC-type star clusters (ONC A in Table[T). Using 
the same assumptions, is lower for all other models in Table [TJ 
due to their higher velocity dispersions. In the case of the UCD- 
models from this paper, the difference between the actual cross sec- 
tion and the geometric cross section is less than 1 per cent. Thus, 
gravitational focussing of stars onto the proto-star plays a minor 
role for the models in TableQ] which justifies the approximation. 

The values for t cnc resulting from equation l[9j are noted in Ta- 
bleland plotted in Figs.Q]and[2] Comparing these values with the 
characteristic life-time of a proto-star, t ploto ~ 10 5 years, it can be 
seen that t cnc < t ploto for the UCD-models with iW*,o = 10 8 Mq . 
For the UCD-models with M*,q = 10 7 Mq and _R p i.o = 3 pc, t onc 
is only slightly larger than t pro to- However, for the UCD-models 
with M*,o = 10 7 Mq and R p i : o = 5pc, t cnc exceeds t pro to by 
about an order of magnitude. This suggests that the encounters be- 
tween proto-stars and stars would influence star-formation in the 
models with M*,rj = 10 s Mq and also, to a much lesser extent, 
in the more compact UCD-models with M*,o = 10 7 Mq, but not 
in the UCD-models with M»,o = 10 7 M and J? p i, = 5pc. The 
UCD-models with M*,o = 10 7 Mq and i? p i,o = 5pc are in this 
respect similar to models A and B for th e initial states of ONC-type 
star clusters from IKroupa et al] d200ll) (see Fig. [2)- This implies, 
invoking the universality of the IMF in open star clusters, that the 
IMF in those UCD-models should also be given by equation 
if deviations from the canonical IMF are caused by encounters be- 
tween proto-stars and stars. According to the calculations in this 
paper, the UCD-models with M».o = 10 7 Mq and fl p i,o = 5pc 
would indeed only evolve into objects similar to an observed UCD 
if their mass loss is as implied by the canonical IMF. The IMF 
in the other UCD-models would however have to be top-heavy to 
some extent for this (see SectionfJJ- The UCD-models in this paper 
are thus self-consistent in that sense. 

We note that equation ^ reveals the particular importance of 
encounters between proto-stars and stars. For a collision between 
two stars, r prot o has to be substituted by a value <C 1 AU, which 
leads to t cnc 3> 10 5 years. Thus, collisions be tween stars only as a 
mechanism that chang es the shape of the IMF (Bonn eiTet al J 19981 
lBonnell&Ba3 ;2002) requires even higher densities. For the en- 
counter between two proto-stars, the density of stars at a given time 
has to be substituted by the density of proto-stars at that time. Tak- 
ing 1 Myr as the characteristic time-scale on which star-formation 
takes place and 10 5 years as the life-time of a proto-star suggests 
that the density of proto-stars is ~ O.lppi,*^, which is five times 
less than the density of stars at the time when half of the total stellar 
population of the UCD has formed. 

A caveat to the above discussion is that it is not specified what 
the consequence of a collision between a proto-star and a star is. 
This is a merger if the encounter is slow enough. If the encounter is 
fast enough for the star to only pass through the proto-star, the star 
transfers some of its kinetic energy on the proto-star and thereby 
disperses some of the matter that would otherwise accrete on the 
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Table 2. The IMFs considered for UCDs. The content of columns is the 
following: Column 1: the identification number of the IMF, Column 2: the 
slope of the high-mass end of the IMF, a (where 2.3 is the Salpeter slope), 
Column 3: the upper stellar mass limit, m max , Column 4: the ratio between 
the total inital mass of stars more massive than 8 Mq and the total initial 
mass of all stars, Column 5: the ratio between the inital number of stars 
more massive than 8 Mq and the initial number of all stars, Column 6: the 
initial mean mass of stars. 



IMF 


a 


'"■max 


A4 mSi o/M, l0 


AT hms>0 /iV„,o 


m 






[M ] 






[M ] 


1 


1.1 


150 


0.921 


0.2031 


10.02 


2 


1.1 


100 


0.886 


0.1830 


7.16 


3 


1.5 


100 


0.719 


0.0632 


2.49 


4 


1.9 


100 


0.453 


0.0210 


1.07 


5 


2.3 


150 


0.230 


0.0072 


0.65 


6 


2.3 


100 


0.213 


0.0071 


0.64 



proto-star. For deciding which of these processes would dominate 
for a given velocity dispersion, as well as for answering the ques- 
tion of how and to what extent they would alter the IMF, detailed 
modelling of the collisions would be required. However, the dis- 
cussion here implies that any process resulting from an encounter 
between stars and proto-stars should only be relevant for the denser 
UCD-models in Tab.[T] in contrast to the models for ONC-type star 
clusters, where most proto-stars should be unaffected by encoun- 
ters. 

We revisit the matter of a possible influence of encounters on 
the IMF in UCDs in Section[3H 



2.1.3 The IMF of the UCD-models 

For each of the nine sets of models listed in Table [T] six IMFs are 
considered. They are either canonical or top-heavy to a different 
degree and have upper mass limits, m max , of either 100 Mq or 
150 Mq. However, all of them agree with the canonical IMF (equa- 
tion QJ for rn < 1 Mq . Studies on m max suggest that m max = 
150 Mq is more realistic than m max = 100 M q for very massive 
star clusters and therefore also f or UCDs (e.g. Mass ev & Hunterl 
ll998tlFiger et alJl998tlFiged2004h . However, the treatment of stel- 
lar evolution and its effect on the mass loss from UCDs in this 
paper (see Section [272b is based on stellar models that only range 
up to a 120 Mq star. Assuming m max = 150 Mq for our models 
therefore requires extrapolating from the given data, which may be 
problematic due to the strong dependencies of stellar properties on 
stellar mass. The emphasis in this paper is therefore on IMFs with 
m max = 100 Mq (which is also a common choice in simple stellar 
population models). The impact of the higher m max = 150 Mq is 
only tested for the canonical IMF and the most top-heavy IMF. 



2.2 Generating the mass loss Tables 

The interstellar medium (ISM) of a new-born star cluster or UCD is 
massively heated by the radiation from massive stars, which leads 
to a mass loss from it until the ISM is depleted. The eventual evo- 
lution of the massive stars into supernovae (SNe) heats the ISM as 
well, but also replenishes the ISM. The rate at which mass is lost 
from the star cluster or UCD due to this interaction between the 
massive stars and the ISM is the driving force for its early evolu- 
tion. This is why the mass loss rate has to be quantified for our 
models. It is recorded in look-up tables, listing how much the mass 
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Figure 3. The data from lSchaller et alj Jl992h on the lifetimes of stars with 
different initial masses from 7 to 120 Mq (open circles) and an interpo- 
lation function to them (solid line), which is given by equation 1121 . It is 
apparent that the stars with the highest masses evolve over an extremely 
short time span. This increases the significance of the upper mass limit of 
the IMF for the dynamical evolution of a star cluster or UCD. 
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Figure 4. The SN rates with time for the modelled UCDs with initial total 
stellar mass M*,o = 1O 7 M . The different curves are for the different 
IMFs listed in Tab. [2] The numbers refer to the labels given to the IMFs in 
Tab. [2] Note that the choice of the upper mass limit of the IMF determines 
the time when the first SN explodes, but turns out to be almost irrelevant 
for the SN rates. For the models with M t>0 = 1O 8 M , the SN rates are 
higher by a factor of 10. The SN rates are propotional to the energy input by 
the SN, because all SN are assumed to release the same amount of energy 
(10 51 erg). 
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Figure 5. The evolution of the stellar mass of the UCD, Mr , with time due 
to stellar evolution. M, is plotted in units of the initial stellar mass of the 
UCD, M„.o- The different lines in this figure represent M* for the different 
IMFs listed in Table E] The high -mass slope of the IMF is noted above the 
curves. Solid (red) lines are for IMFs with an upper mass limit of 100 Mq 
and dashed (blue) lines are for an upper mass limit of 150 Mq. The choice 
of the upper mass limit has only a minor impact for the canonical IMF 
(a = 2.3), but is more significant for the most top-heavy IMF we consider 
(a = 1.1). For models with SFE=1 and HE=1 these curves show the total 
mass loss as well. 



of the UCD-models has to be reduced for each time-step in the cal- 
culation. 

Evidently, knowing the lifetimes of massive stars is essen- 
tial for generating the mass loss tables. A very good proxy for the 
time at which the life of the star ends is the time at which carbon 
burning has finished. The time that has elapsed until this evolution- 
ary stage is reached is taken from the stellar evolutionary grid by 
Schall er et all i 19921) for massive stars with various initial masses. 
This time-span is identified with the lifetime of a star in this paper. 
A good fit to the lifetimes of stars with high mass (m, ^ 7M0) 
and low metallicity (Z — 0.001 and [Z/H] = —1.3 respectively) 
is the function 

m, = a(t, - b) c , (12) 
with 

a = 74.6, 6 = 2.59, c = -0.63. 

where the initial mass of the star, m,, is measured in Mq and the 
lifetime of the star, t,, is measured in Myr (Fig. [3j- It thus cov- 
ers the whole range of stellar initial masses of stars that undergo 
SN explosions at the end o f their evolution, which is m, > 8 Mq 
dKoester & Reimers|[l99r3) . The parameters a, b and c have been 
found by a least-squares fit. The models for low-metallicity stars 
were preferred over models for stars with Solar metal licity because 
of the mostly sub-solar meta llicities of the UCDs I IMieske et alj 
l200dlEvstigneeva et al.l2fJ07h . This choice has however only a mi- 
nor impact on the best-fitting parameters a, b and c. Note that since 
all stars in the UCD-models are assumed to have formed at once 
in this paper, it is possible to substitute the stellar lifetime, t,, in 
equation U2\ with the age of the UCD-model, t, in order to find 
the initial mass of the stars that undergo SNe at that time. 

Now consider the increase of the age of the UCD-model by 



the time step ti ti+i. During this time, AiV» stars with a total 
mass AM* will complete their evolution. These quantities can be 
written as 



AN*a 



and 



AM, i 



M» 



Mr. 



M*. 

~Wr. 



£(m*) dm, 



^(m,)m, dm,, 



(13) 



(14) 



where £(m») is the IMF, m*,i is the initial mass of stars that evolve 
at t = ti and m*,i+i is the initial mass of stars that evolve at t = 
ti+i. M,fi is the total initial stellar mass of the UCD-model. Given 
the normalisation chosen for the IMF (see equations Q] and |2), the 
purpose of the factors M*,o /Mq is to scale equations |13IU4l and|15l 
to a UCD-model with the initial mass of M* j0 . AN, is equivalent 
to the number of SNe during the time step i; — > ti+i; i.e. the SN- 
rate in the limit of ti+i — ti — ► 0. At the time when the most 
massive stars evolve, this SN-rate is, for instance, ~ 1 SN per 10 
years for the UCD-models with iW*,o = 10 7 Mq and a high-mass 
IMF slope of a — 1.1, while it is a few SN per 10 3 years for the 
UCD-models with M*,o = 10 7 M and a = 2.3. The influence 
of the top-heavyness of the IMF on the SN-rates decreases as time 
proceeds. The SN-rates for the UCD-models with M*, = 10 7 M 
and the IMFs from Table [2] are shown in Fig. [4] The SN-rates for 
the UCD-models with M, fl = 10 s M Q are higher by a factor of 10 
compared to the ones shown in this figure, but the same otherwise. 

Fig.[5]depicts the change of the stellar mass of the UCD-model 
with time, i.e. M*,t = A4»,o for ti ^ *sn and M*,s = M*,o — 
Yy„- 1 AM*, n for ti > <sn, where M*,o is the total initial stellar 
mass and tsN is the time when the first stars become SNe. 

The total energy deposited by stars into their surroundings by 
radiation and stellar winds at the time t = ti, -L*,i, is given as 



L,s = 



M, 



,0 



£(m*)Z(m«) dm,, 



(15) 



with m max ,i being the mass of the most massive star that has not 
evolved into a SN at that time. Z(m*) is the energy deposition rate 
of stars into the ISM through radiation and stellar winds as a func- 
tion of their initial mass. It is estimated as 



Z(m») = 2.16 x 10 



(16) 



m, \ erg 
Mq"J Myr' 

w hich is identical to equa tion (12) inlBaumgard t et al.l d2008l) . As 
in lBaumgardt et ail d2008l) . equation J 1 6b is applied to stars of all 
masses, even though it was obtained in a fit to high-mass stars. 
Note that the positive exponent in equation J16l > and the negative 
exponent in the IMF cancel out more or less in equation dl5t . The 
contribution of low- and intermediate-mass stars to the total en- 
ergy deposition into the ISM is therefore small at first, because the 
masses of high-mass stars are distributed over a much wider range. 

2.2.7 The algorithm 

The integrations in equations {Yi\ to d 1 5 are done numerically for 
the IMFs listed in Table l[2j. The used program is structured as out- 
lined below. 

Start at t = with a set of initial parameters taken from Ta- 
ble[T]and an IMF taken from Table|2] Let At be the time step from 
ti to ti+i and tsN the time when the first stars become SNe. 

(i) If ti > tsN, calculate which stars evolve from t = ti to t = 
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Figure 6. The change of the total mass of the UCD-model, M p i, with time due to stellar radiation and evolution for our models with SFE = 0.4 and 
HE = 1.0 or HE = 0.03 relative to the total initial mass of the UCD-model, M p i t Q. The assumptions regarding SFE, HE, initial mass of the UCD-model and 
inital Plummer-radius of the UCD-model for the mass loss histories shown are indicated in the corresponding panel. The different lines in this figure represent 
M p i for the different IMFs listed in Table E] The mass loss increases with the top-heavyness of the assumed IMF: for the topmost curves a = 2.3 and for 
the lowermost curves a = 1.1. Solid (red) lines are for IMFs with an upper mass limit of 100 Mq and dashed (blue) lines are for an upper mass limit of 
150 Mq. The choice of the upper mass limit has only a minor impact for the canonical IMF (a = 2.3), but is much more significant for the most top-heavy 
IMF that are considered (a = 1.1). Note that for t > 3 Myr the mass loss histories shown in this figure are equal to the change in stellar mass shown in Fig. [3] 
if the primordial gas is expelled before the first star has evolved completely. However, contrary to Fig. |5]the time-axis is scaled logarithmically here, in order 
to show the sometimes very rapid expulsion of the primordial gas. 



ti + At using equation dl2t and then which total mass these stars 
have, AM„,i (equation [14} • This mass is added to the total mass of 
the interstellar medium, MisM.i- This includes the possibility that 
the UCD-model had no ISM left at the end of the previous time 
step. In this case, MisM.i = AM,j. 

(ii) The rate at which the stars and the SNe deposit energy into 
the ISM during the time step is calculated, Li. Li is approximated 
by Li = L»,i+i + LsN,i, where L„,i+i is the rate at which the 



stars deposit energy into the ISM at t = ti+i (equation 1 1 5 1> and 
LsN.i is the energy that the SNe deposit into the ISM from t — U 
to t — ti + At (which is in the limit of At — > an energy de- 
position rate as well). Using the number of stars that evolve dur- 
ing the time step AiV,,!, (equation [T5), LsN,i can be estimated by 
assuming that each SN releases a characteristic amount of kinetic 
and electromagnetic energy, which are the forms of energy that are 
relevant for driving matter out of the UCD-model. Estimating this 
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quantity as 10 51 erg per SN (e.g. ICarroll & Ostlid[l996h leads to 
Lsn,i = 10 51 AN, t i erg Myr . The total luminosity is multiplied 
by the HE to obtain L]d n ,i, which is the luminosity that is not ra- 
diated away through thermal emission of the ISM, but is converted 
into kinetic energy of the gas leaving the UCD. 

(iii) The time tj it would take until all gas is expelled from the 
UCD-model is estimated, assuming that Lkm.i does not change dur- 
ing that time. This is done using the equation Til/kin, i = \E po t t i — 
E po t*,i\, where -Epot.i is the total binding energy of the UCD- 
model at t = U and £/ po t*,i is the binding energy the UCD would 
have if it would lose all gas at that time. Note that the UCD-model 
inflates as it loses mass and -B po t*,i should therefore be calculated 
using the Plummer-radius the UCD-model has after all gas is ex- 
pelled. We estimate it using the relation between initial radius and 
final radius of a stellar system for adiabatic mass loss, even though 
the mass loss is in our case sometimes clearly not adiabatic. This 
relation is given by 

rfinal _ Mjnit „_ 
Pinit A/final ' 

where rjrjt and Afjnit are radius and mass of the stellar system at 
the beginning of mass loss, respectively, and rfl na i and Mfi na ] the 
are radius an d mass of the stellar system at the end of mass loss, 
respectively (Kroupa 2008 and references therein). Equation dl7t 
underestimates the expansion of the UC D-models at ti mes of non- 
adiabatic mass loss (cf. equation 8.20 in lKroupa|[2008l) . Therefore, 
if the UCDs have experienced extended non-adiabatic mass loss, 
the mass loss rates calculated here are too low, because a more pro- 
nounced expansion implies that the potential well becomes shal- 
lower and the remaining gas requires less energy to escape from it. 
This is even more important for a more top-heavy IMF. 

(iv) If n < At, set Ti = At. The mass loss of the UCD-model 
during the time step is assumed to be 

SM = M lS M,i — ■ (18) 

Ti 

(v) Calculate the new parameters of the UCD-model after it has 
lost the mass SM : The new total stellar mass is decreased by AM*; 
and the new Plummer-radius is estimated using equation jl7\ . 

(vi) If ti+i = ti + At is less than it takes a star with m* = 
8 Mq to evolve into a SN according to equation d!2t . repeat steps 
(i) through (vi), but for t;-|_i instead of ti. 

The underlying assumption in the chosen approach is that the 
material expelled from a SN does not immediately escape the UCD, 
but that its kinetic energ y is thermalised, as i t is a ssumed for mas- 
sive star clusters in, e.g., lTenorio-Tagle et al.l J2007h . This can hap- 
pen either through interaction with the surrounding ISM or through 
the collision of the expanding envelopes of different SNe. The latter 
becomes more relevant with increasing top-heaviness of the IMF. 
The notion of the thermalisation of the SN ejecta is flawed if there 
is no ISM left and if the SNe are too few for their envelopes to in- 
teract with one another. However, in this case also very low HEs 
are sufficient to keep the UCD-models gas-free. 

The mass loss histories calculated by using the above routine 
are shown in Fig. ((6} for the models with SFE=0.4 and HE=1 or 
HE=0.03. For the the models with SFE=1 and HE=1, the evolution 
of the stellar mass of the UCD-models shown in Fig. [5] also illus- 
trates their mass loss history, since the UCD-models are gas-free at 
all times in this case. 



Table 3. The total masses of all stars more massive than 8 Mq (Mh ms o, 
Column 2) and the total masses of their compact remnants (M rem , 
Columns 3 and 4) for the IMFs in Tablef2] The masses are in units of the to- 
tal mass of all stars that were formed initially, M*.o - The two different val- 
ues for Mrem for a given IMF reflect that the mass of the remnants of very 
massive stars is poorly known. While the mass of the compact remnants of 
stars with initial masses of 8Mq ^ m, < 25 Mq is 1.35 Mq in both 
estimates, the mass of the compact remnants of stars with m* > 25 Mq is 
assumed to be either 0.1m* (Column 3) or 0.5m* (Column 4). 



IMF 


A/hms,0/jW*,0 


M rcm /M*, 
tobh = 0.1m* 


Mrem/M*,o 
tobh = 0.5m, 
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0.921 


0.0910 


0.409 
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0.0871 


0.369 
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0.719 


0.0709 


0.271 
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0.453 


0.0452 


0.150 
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0.230 


0.0234 


0.0698 
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0.0218 


0.0605 



2.2.2 The role of compact stellar remnants 

A simplification that is made in the creation of the mass loss histo- 
ries is that the whole mass of the evolved stars is added to the ISM, 
including the mass of their compact remnants. For testing under 
which conditions this approximation is reasonable, the total mass 
of all compact remnants of stars with initial masses m* > 8Mq, 
M rem , needs to be compared to the total mass of their progenitors, 
Mhms,o- If an IMF is given, calculating Mrem requires a relation 
between the initial masses of stars and the masses of their compact 
remnants. Such an initial-to-final mass relation is, e.g., formulated 
in equation (8) of DKB. Their equation is also used here. Thus, 
stars with initial masses of 8 Mq ^ m* < 25 Mq are thought 
to evolve i nto neutron stars (NSs) with a m ass of 1.35 Mq, which 
is the mass iThorsett & C hakrabartvl dl999h have found for pulsars, 
i.e. a sample of neutron stars that can easily be detected. Stars with 
even higher initial masses are believed to evolve into black holes, 
but the actual masses of thes e black holes (BHs) are poorly con- 
strained (cf. figs. 12 and 16 in lWooslev et al.l2002h . Therefore, two 
cases are considered for the masses of BHs, namely the case that 
they all have 10 per cent of the initial mass of their progenitors 
(mBH = 0.1m*) and the case that they have 50 per cent of the 
initial mass of their progenitors (niBH = 0.5m,) are considered. 
The resulting values are noted in Table [3] 

It is apparent from these numbers that for mBH = 0.1m, the 
mass locked up in compact remnants is indeed negligible, while 
this is not the case for mBH = 0.5m*. Howev er, the masses 
of observationally confirmed stellar-mass BHs (see lCasaresll2007l) 
seem to favour the case of mBH = 0.1m*, leadi ng to BH masses 
< 10 Mq. Apart from that. lLvne & LorimeJ Jl994f) report a mean 
birth velocity of 450 ± 90kms _1 for pulsars (i.e. neutron stars) 
and the processes that precede the birth of a stellar mass BH are 
esse ntially the same as th e ones that precede the birth of a neutron 
star dWooslev et aiT2 002). This suggests that a large fraction of the 
compact remnants (BHs as well as neutron stars) are born with ve- 
locities well above the escape velocity of the UCD-models in Tab. 1, 
which for a Plummer sphere is about twice the ve locity dispersion 
(compare equations 8.59 and 8.61 in lKroup3l2008l) . Thus, in a re- 
alistic scenario, the total mass of the compact remnants remaining 
in the UCD is likely to be small compared to the total mass of the 
progenitors of all compact remnants. Moreover, the mass-loss his- 
tories created by the algorithm described in Section |2.2. 1 1 suggest 
that the UCDs are gas-free at the end of the evolution of massive 
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stars, with the exception of the models with high initial mass and 
low heating efficiency. The latter models suggest however that the 
UCDs consist mainly of gas at that time, which seems unlikely (see 
Section [3~3l >. As a conclusion, UCDs are likely to have lost most of 
the mass that was locked in massive stars at the time when massive 
stars have evolved, if they formed as is assumed here (i.e. as very 
massive star clusters). This mass loss proceeds however not only by 
the escape of the gaseous components of the SN-remnants from the 
UCD (i.e. a process modelled by the algorithm in Section |2. 2. U , 
but also by the ejection of the compact remnants. This latter pro- 
cess would play a substantial role for the total mass loss of the UCD 
if niBH = 0.5m*, but not if 771bh = O.Ira*. 

Note that the expectation of a large difference between the 
total mass of the compact remnants left in UCDs and the total mass 
of their progenitors is also the reason why the IMFs of the UCDs 
have to be so extremely top-heavy, if the enhanced M/Ly ratios 
of the UCDs are to be explained by an over-abundance of stellar 
remnants. This is what motivated the sometimes extreme choices 
for the IMF in the UCD-models in the first place. 

In essence, neglecting the remnant masses seems justifiable 
in the context of the present study, as it also helps to avoid a num- 
ber of very speculative assumptions. This includes the precise mass 
of the remnants, which fraction of them remains in the UCDs and 
how much of the kinetic energy available from the SNe is tran- 
ferred to them. As a result, the mass loss is over-estimated in the 
UCD-models, but probably not by much more than 10 per cent. 
Consequently, their expansion is over-estimated as well by about 
the same amount, if the heating is sufficient to expel all gas from 
them. This bias is thus opposed to the bias induced by the assump- 
tion of adiabatic mass loss at all times in the calculation of the mass 
loss histories. 

We note that also the treatment of the energy input from SNe 
(each of them contributes 10 51 erg) and the energy input from stars 
(equation[T5]using equation[T6t is only approximate, but can hardly 
be done with greater precision with current knowledge. 

2.3 Time evolution of the UCDs 

The UCD-models are set up to be in virial equilibrium before the 
onset of mass loss. This is motivated by the fact that star-formation 
in a star cluster takes place on a time-scale of w 1 Myr, while the 
crossing times in the UCD-models are about an order of magni- 
tude lower (see Table [T) a nd the time-scale for v iolent relaxation 
is a few crossing times i Bi nnev & Tremainell 19871) . Thus, the time- 
scale for the UCD-model to settle into a state near virial equilibrium 
is shorter than the time-scale for the formation of its stellar popu- 
lation. Note that the assumption of virial equilibrium is crucial for 
the validity of the results in this paper, since UCDs would evolve 
completely different if the y were no t in vi rial equilibrium at the 
onset of gas expulsion, see iGoodwinl 020091 ). But it is also argued 
there that very massive star clusters are much more likely to be in 
virial equilibrium at that time. 

To calculate the evolution of the UCDs in the fir st few Myr the 
particle-mesh code Superbox dFellhauer etafl boOO) is used. Each 
UCD is represented by 1 million particles and is integrated forward 
in time until 200 Myr using a small time-step of 0.01 Myr for the 
models with HE=1 and SFE=1 and a time-step of 0.005 Myr for all 
models with SFE=0.4. The smaller time-step for the models with 
SFE=0.4 is necessary because of their shorter crossing times due 
to their higher initial masses for our assumed stellar masses, see 
Table Q] The code is altered to allow for the mass loss due to gas 
expulsion and rapid stellar evolution in the first tens of Myr. To 
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Figure 7. Change of the Lagrangian radii (10, 20, ... 90 % mass) with time 
for model m8j-5.sl.hl with IMF 1 (m ma x = 150 Mq, a = 1.1). 

mimick this mass loss we implemented the look-up tables whose 
generation is described in Section I2T2I They give the total mass of 
the UCD at each time-step. The mass of each particle and hence- 
forth the total mass of the modelled UCD is reduced accordingly. 

The UCDs are modelled in isolation, i.e. in the absence of a 
tidal field, even though UCDs are found in the vicinity of mas- 
sive elliptical galaxies. But regarding the short time span of our 
computations of 200 Myr (compared to their orbital times of a 
Gyr or longer) and the fact that e.g. at a distance of 80 kpc and 
adopting the potential of M 82 the tidal radii would be 600 pc 
for the models with M p i,o = 10 7 M Q and 1400 pc for the mod- 
els with Afpi^o = 10 s Mq, th e effect of the tida l fields can be 
neglected. (See also tab le 6 in iHilker et al.l 120071 and table 8 in 
lEvstigneeva et alj|2007l for a comparison between half-light radii 
and tidal radii of UCDs.) 



3 RESULTS 

We calculated a suite of 56 models, combining each of the sets of 
UCD parameters given in TableQ]wifh each IMF given in Table|2] 
The results are discussed separately for the different assumptions 
regarding the SFE and the HE in the following, see also Tables 3, 4 
and 5. 

3.1 SFE=1 

The assumptions SFE=1 and HE=1 stand for the case of highly 
efficient star-formation and heating. There is no expulsion of pri- 
mordial gas in this case, but the mass loss through the evolution of 
massive stars can still be quite severe (as the UCD is cleared eas- 
ily from the products of stellar evolution with such a high HE). It 
amounts to up to about 90 per cent of the initial mass for the UCDs 
with the most top-heavy IMFs. However, this mass loss is slow 
compared to the short crossing-times of the initially very massive 
and compact models. This makes sure that the calculated UCDs al- 
ways survive this period of mass loss. In Fig. [7] the time evolution 
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Table 4. Final results of the calculations for a SFE of 1 and a HE of 1 . The models whose final parameters match the observed parameters of UCDs best 
are marked with a (+) before the first column. The information given in the columns is the following: Column 1: the name of the model as given in Table[T] 
Column 2: the IMF as given in Table|2] Column 3: the total mass of stars that have not evolved at the end of massive-star evolution (i.e. stars with m < 8 Mq) 
in units of the initial total mass of all stars, Column 4: the total mass of stars that remain bound to the cluster in units of the total mass of all stars less massive 
than 8 Mq (i.e. no stars become unbound if the entry in this column is 1), Column 5: the mass of the cluster at the end of the calculation, Column 6: the 
final half-mass radius, Columns 7 and 8: the final Plummer-radius fl p i,f and its 1-cr error, Column 9: the expansion factor / c , Columns 10 and 11: the central 
surface density £ f with its error, Columns 12 and 13: and the central line-of-sight velocity dispersion cr f with its error. E £ and ctq £ are derived by fitting 
Plummer-profiles to the data at t = 200 Myr, using a non-linear least-squares Marquardt-Levenberg algorithm. Both the fits to £q f and to cfq f also deliver 
estimates for iJ p i,f . The quoted value for Rpif is the one obtained from the fit to £o,f . but the one obtained from the fit to crrj,f is not much different. 
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Figure 8. Expansion factors f c of our models for a star formation effi- 
ciency of 1 and a heating efficiency of 1 plotted against the number given 
to the assumed IMF (as in Table[2}. The symbols show the different initial 
structural parameters of the UCD-models: (red) crosses for Rpifi = 3pc 
and Mpi o = 10 7 Mq, (green) squares for J? p i,o = 5pc and Af p i o = 
10 7 Mq and (blue) circles for R plfi = 5 pc and M ply0 = I0 7 M Q . 



of the Lagrangian radii of one of the models is shown. It can be seen 
clearly that after an interval of rapid expansion due to the mass loss 
the UCD-model finally settles back into a new equilibrium. The ex- 
pansion factor / e of the models is measured by comparing the final 
(-Rpi.f) with the initial (i? p i,o) Plummer-radius, 



fc = 



-Rpi.o 



(19) 



The Plummer-radii are found by fitting Plummer-models to the sur- 
face density profiles of the UCDs, using a non-linear least-squares 
Marquardt-Levenberg algorithm. The Plummer-radius is also iden- 
tical to the projected half -light radius for the models (see equa- 
tion 8.57 in lKroupall2008l) . Fig. [8] shows the expansion factors for 
all UCD-models with SFE=1 and HE=1. It is visible that among the 
clusters with top-heavy IMFs (IMFs 1 to 4) the UCD-models with 
the highest mass expand the least. This is because the more massive 
UCD-models have shorter crossing times and are therefore closer 
to the regime of adiabatic mass loss. 

Table|4]shows the final quantities for the models with SFE=1 
and HE=1. The models that are the best representations of present- 
day UCDs at the end of the calculation are marked with a '(+)' in 
front of the first column. Note that some, but not all of these have 
the canonical high-mass IMF slope. 

We note that assuming HE=1 is not decisive for most of the 
UCD-models with SFE=1. For instance, if HE=0.03 is assumed for 
the UCD-models with an initial stellar mass, M,, ,of 10 7 M Q and 
SFE=1, the mass-loss histories of such models are the same as in 
the case of HE=1. The same is true for UCD-models with 10 8 Mq 
SFE=1 and HE=0.03 if their IMF is the canonical one. On the other 
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Figure 9. Expansion factors / c of the models with a star formation effi- 
ciency of 0.4 and a heating efficiency of 1 plotted against the number as- 
signed to the assumed IMF (as in TablefS}. The symbols show the different 
inital structural parameters of the UCD-models: (red) crosses for R p i t o = 
3pc and M p \ q = 2.5 X 10 7 Mq, (green) squares for R p ifi = 5pc 
and Mpi o = 2.5 X 10 7 Mq and (blue) circles for R p i t o = 5pc and 
Af plj0 = 2.5 X 10 7 M . The modelled UCDs with IMFs 1 and 2 dissolve 
completely due to their heavy mass loss. 

hand, for the two most top-heavy IMFs in Tab. [2] (IMFs 1 and 2, 
q = 1.1), the UCD-models with 10 8 M and SFE=1 retain most 
of the gas released by the evolution of massive if HE=0.03, while 
they are gas-free at all times if HE=1. In the case of H=0.03, these 
models are very similar to the UCD-models with 10 8 Mq, SFE=0.4 
and HE=0.03, which are discussed in Section [3~3l 

Thus, it is UCD-models with the most top-heavy IMFs that re- 
tain gas the easiest. This is because by assuming that the amount of 
energy released by a SN does not depend on the mass of the progen- 
itor star (as done in this paper), the total mass set free by the SNe 
increases more quickly than the total energy provided by the SNe. 
Also, the luminosity of the stellar population of the UCD, which is 
the other energy source that powers its mass loss, decreases more 
rapidly with time for more top-heavy IMFs. 

3.2 SFE=0.4 and HE=1 

The assumptions SFE=0.4 and HE=1 imply an even more dramatic 
mass loss than the case of SFE=1 and HE=1 (Section r3.lt . The en- 
ergy input of massive stars is high enough to clear the UCD-models 
of the primordial gas either well before or at the time the first stars 
end their evolution on the main-sequence. The mass loss is in fact 
so rapid that a significant fraction of the stars of the UCD-models 
become unbound, even if an IMF with the canonical high-mass 
slope (q = 2.3, IMFs 5 and 6) is assumed. The calculated UCD- 
models dissolve completely if the IMFs with the flattest high-mass 
IMF-slopes (a = 1.1, IMFs 1 and 2) are assumed. Note that they 
are dissolved by the combination of the very rapid expulsion of the 
primordial gas and the more gentle mass loss though stellar evolu- 
tion, since an instantaneous loss of 60% of the init ial mass would 
still leave a bound remnant dBoilv & Kro"upall2003l) . as would the 
mass loss through stellar evolution alone (cf. Section [3~ll . 




t [Myr] 

Figure 10. As Fig. El but for model m8_r5.s04_hl with IMF 4 (m max = 
100 M Q , a = 1.9). 

Analogous to Section [3~T1 the expansion of the UCD-models 
is measured by the ratio between their final Plummer-radii and their 
initial Plummer-radii and the results are plotted in Fig. [9] It turns 
out that in this set of models, the UCDs that expand the most are 
always the ones with the longest crossing times while the UCDs 
that expand the least are always the ones with the highest initial 
mass. 

Table[5]shows the final quantities for the models with SFE=0.4 
and HE=1. Model m8_r5_s04.hl with IMF 4 (a = 1.9) is the only 
one of them with a top-heavy IMF and with a good agreement be- 
tween its final parameters and the parameters observed in UCDs. 
The evolution of its Lagrange-radii is shown in Fig. 1101 

3.3 SFE=0.4 and HE=0.03 

In the case of a moderately high star formation efficiency 
(SFE=0.4) and and low heating efficiency (HE=0.03), the UCD- 
models with Mpi,o = 2.5 x 10 7 Mq are gas-free at the end of the 
computation (as are all models with HE=1). In contrast to that, the 
models with M p i,o = 2.5 x 10 8 Mq keep most of the gas at such a 
low HE, so that these UCD-models are predicted to expand barely 
and to consist mainly of gas at the end of the integration, implying 
very high M/Ly ratios. (In the two most extreme cases, where the 
high-mass IMF slope is a = 1.1, approximately 5% of the total 
mass of the UCD-model is stars at that time, while the rest is gas.) 
This is a very implausible situation. It is more likely that if the heat- 
ing efficiency is too low to drive the gas out of the cluster, the star 
formation efficiency would become higher through new star forma- 
tion episodes, until eventually all matter is locked up in low-mass 
stars and thereby a SFE of 1 is approached. 

However, the half-mass radius of a UCD would hardly change 
with time in this case, since it keeps most of its initial mass (cf. 
equation [TTt. while the initial half-mass radii of UCDs suggested 
in this paper are clearly smaller than t he half-mass radii of ob- 
served present-day UCDs (cf. table 5 in lMieske et alj l2008). This 
means that, if a UCD indeed retains most of its mass, it must be 
born with an initial half-mass radius close to the observed values. 
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Table 5. As Table [4] but for a SFE of 0.4 and a HE of 1. These models predict the complete dissolution of the UCD if IMF 1 or IMF 2 are assumed (the ones 
with high-mass slope a = 1.1). 
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Table 6. As table |4] but for a SFE of 0.4 and a HE of 0.03. These models predict the complete dissolution of the UCD if IMF 1 or IMF 2 and an initial total 
mass of 2.5 X 10 7 are assumed. The models with an initial total mass of 2.5 X 10 s Mq on the other hand keep most of their primordial gas and therefore 
their evolution is completely different. 
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By calculating mass-loss histories of UCD-models with M p i )0 = 
2.5 x 10 s M Q , SFE=0.4 and HE=0.03 for different initial Plum- 
met' radii (by the method described in Section [2.2. lb . it turns out 
that the energy input from massive stars is sufficient to remove all 
gas from these UCD-models at an initial projected half-mass ra- 
dius (i.e. Plummer radius) of 12 pc instead of 5 pc, even if the 
UCD-models have the canonical IMF. In contrast to that, the most 



massive observed UCDs, with masses of « 10 Mq, are reported 
to have half-mass radii of ~ 100 pc. Thus, adopting a major star 
burst as the scenario for the birth of a UCD (as done in this pa- 
per), an object that is able to keep gas after the star burst would be 
too compact to evolve into an UCD. On the other hand, if the object 
has an initial half-mass radius that allows it to evolve into a UCD, it 
would loose its gas on a time scale of a few Myr. This excludes the 
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Figure 11. As Fig. [9] but for a heating efficiency of 0.03 instead of 1. The 
modelled UCDs with M» i0 = 10 7 M and IMFs 1 or 2 dissolve com- 
pletely due to their heavy mass loss and are therefore not shown here. 
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Figure 13. Change of the Lagrangian radii (10, 20, ... 90 % mass) with time 
for model m7.r3.s04_h003 with IMF 1 (m max = 150 M , a = 1.1). 
Here the UCD is disrupted by its mass loss, unlike the models shown in 
Figures.|7l[T0land[T2l 
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Figure 12. As Fig.|7] but for model m7j-3.s04.h003 with IMF 4 (m ma 
100 M©, « = 1.9). 



above scenario where a UCD forms a substantial part of its stellar 
population over a longer period of time after the initial star burst. 

Fig.[TT]depicts the expansion rates of the UCDs with SFE=0.4 
and HE=0.03 and shows a strong difference between the modelled 
UCDs with M p i,o = 2.5 x 10 8 M Q and the ones with M p i,o = 
2.5 x 10 7 M . While the less massive UCD-models expand almost 
as much as in the case of HE=1 (see Sections |3. 1 1 and [3.2t , the 
extension of the more massive UCD-models hardly changes. 

The final quantities found for the models with SFE=0.4 and 
HE=0.03 are shown in Table [6] Except for models with IMFs that 
have the canonical high-mass slope, only model m7_r3_s04_h003 
with IMF 4 (a = 1.9) is a good representation of a UCD. The 



time-evolution of its Lagrange-radii is shown in Fig. [12] Fig. [13] 
on the other hand illustrates the evolution of a UCD that dissolves 
completely because of extreme mass loss due to its very top-heavy 
IMF. The only difference to the model shown in Fig. [T2] is that 
IMF 1(« = 1.1) instead of IMF 4 was assumed. 



3.4 Implications on the initial parameters of UCDs 

Based on the fraction of the mass that is lost from the modeled stel- 
lar systems and the factors by which they expand due to mass loss, 
initial conditions that would lead to UCD-like objects can be esti- 
mated. The results of such estimates are shown in Fig. [14] Thus, 
UCDs may have been born from extremely compact configurations 
with densities ranging up to lO 8 M0pc~ 3 . These numbers have 
admittedly to be taken with caution, since the expansion factors 
and the total mass loss of the objects have been derived for stellar 
systems with different initial parameters, using mass loss histories 
through stellar processes that were specifically created for them. 
Note however the similarity between the expansion factors of mod- 
els with the same initial mass and IMF, but different initial radii 
(Figs. [8ll9l and fTTTl. Analogous calculations to the ones performed 
here, but with the initial parameters plotted in Fig.[l4]are therefore 
likely to lead to final parameters that represent the actual parame- 
ters of UCDs better, but this needs to be studied in follow-up work. 

The initial central densities following from the pairs of initial 
masses and initial Plummer-radii plotted in Fig. [14] are shown in 
Fig.[l5] The initial parameters that would lead to UCD-type objects 
according to Figs.[l4]and[T5]can be compared to the initial parame- 
ters of the UCD-models listed in TabQ] whose early evolution was 
calculated in this paper. It thereby becomes apparent that the initial 
conditions resulting in UCDs may be even more extreme than the 
ones that are specified in Table Q] Thus, encounters of proto-stars 
with stars, as discussed in Section [2.1.2| may be even more relevant 
for the star formation in actual UCDs than for the UCD-models cal- 
culated in this paper. The UCDs may even have been dense enough 
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Figure 14. Estimated initial masses and Plummer-radii that would lead to 
representative, UCD-type objects. The black square to the left represents in 
each panel an observed typical small UCD with a mass of 5 X 10 6 M Q 
and a Plummer radius of 10 pc, whereas the black square to the right rep- 
resents an observed typical large UCD with a mass of 40 x 10 6 M Q and 
a Plummer radius of 20 pc. The remaining symbols show estimated initial 
masses and Plummer-radii, that would lead to one of these two represen- 
tative UCDs with the IMFs from Table [2] identified here by the number 
assigned to them in that table. The assumed star formation efficiency and 
heating efficiency are indicated at the top of each panel. For the more mas- 
sive UCD-like object, the estimated initial parameters are based on the total 
mass loss (through stellar and dynamical evolution) and expansion factors 
of the models starting with a total stellar initial mass M„ ,o = 10 s Mq (Ta- 
ble [T}, while the estimates for the less massive UCD-like object are based 
on the models with M„ o = 10 7 Mq (Table[T}. The dashed lines in each 
panel show constant central densities, starting from 10 3 Mq pc — 3 and in- 
creasing by a factor of ten downward with each line. Note that the initial 
conditions resulting in the more massive representative UCD in the lower- 
most panel (SFE=0.4 and HE=0.03) are based on models where hardly any 
gas is lost from the UCD-model while no more stars are formed from this 
material, which is a unrealistic scenario (see Section |33V 
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Figure 15. The initial central densities that would lead to the representative 
UCD-type objects shown as black squares in Fig. [14] given the mass losses 
and expansion factors of the UCD-models calculated in this paper for the 
IMFs listed in Table|2] It is assumed (as in Fig. 1141 that the mass loss and 
the expansion experienced by an object that evolves into the more massive 
representative UCD (with a mass of 40 X 10 6 Mq and a Plummer radius of 
20 pc) is given by the mass losses and expansions calculated for the UCD- 
models that start with an initial total stellar mass of 10 s Mq (Table [T}. 
The evolution of the less massive representative UCD (with a mass of 5 X 
10 6 Mq and a Plummer radius of 10 pc) from its initial state is thought 
to be consistent with the mass losses and expansion factors found for the 
UCD-models starting with an initial total stellar mass of 10 7 Mq (TableQ]. 
The central densities of the initial states of UCDs in this figure are thus 
given by the initial masses and Plummer radii assigned to them in Fig. H4\ , 
using equation (4) with R = 0. 



for frequent collisions bet ween stars, so that thi s process could also 
have shaped their IMF (cf. lBonnell et al.ll 19981) . 

Note the similarity of Figs. ([Sj l[9j and i ll It with the corre- 
sponding panels of Fig. dl5t . except for the different scaling. This 
is because the expansion factor enters with the third power into the 
calculation of the initial density for a given final mass and final 
Plummer-radius according to equation (4j, while the dependency 
on the lost mass is only linear. 
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Also note that the negligence of compact remnants induces 
a bias on the estimated initial parameters: Remnants kept by the 
UCD-model diminish the mass that leaves the UCD-model and 
thereby also reduce its expansion. For arriving at the mass and the 
Plummer radius of the representative UCD-type objects plotted in 
Fig-UU a UCD-model that keeps some of the mass of the massive 
stars in the form of remnants would thus need a larger initial ra- 
dius and a smaller initial mass than a UCD-model that looses all 
remnants from massive stars. Consequently, the initial density of 
the UCD-model that keeps some remnants would also be smaller. 
The total mass of the remnants remaining in the UCDs has however 
been argued unlikely to be much larger than 10 per cent of the to- 
tal mass of their progenitor stars (Section[2j2j2} and the bias on the 
initial parameters shown in Fig. [14] would be of the same order. The 
large implied mass-loss through the evolution of massive stars does 
not contradict the high M/Ly ratios of UCDs, if the number of 
massive stars was sufficient in them, i.e. their IMF was top-heavy 
enough (cf. DKB). 

The bias caused by the the negligence of compact remnants 
may however be alleviated by an opposed bias. This opposed bias 
comes from the fact that non-adiabatic behaviour of the UCD- 
models was taken into account in the actual calculation of their dy- 
namical behaviour, but not in the modelling of the mass loss driving 
the early evolution of the UCD-models (see Section |2.2.1| >. 



4 SUMMARY AND CONCLUSIONS 

We calculate the early evolution of extremely massive star clusters, 
using the particle-mesh code Superbox. Their initial radii are cho- 
sen in concordance with typical values for globular clusters (GCs), 
while their initial masses reflect the masses of ultra compact dwarf 
galaxies (UCDs). The early evolution of a star cluster is driven by 
mass loss through gas expulsion and stellar evolution. This mass 
loss is treated by reducing the mass of each particle in accordance 
with previously tabulated mass loss histories, so that the total mass 
of all particles agrees with the total mass of the UCDs as given 
in those tables. The rate and the magnitude of the mass loss de- 
pends in particular on the stellar initial mass function (IMF). Since 
it was suggested that UCDs may have formed with a top-heavy 
IMF (DKB), the integrations use mass loss tables not only for the 
canonical IMF but also with different top-heavy IMFs. A possible 
explanation for why the IMF in UCDs could be top-heavy is en- 
counters between proto-stars and stars. If UCDs indeed formed as 
the most massive star clusters, as suggested in this paper, such en- 
counters would be quite likely in emerging UCDs. In contrast, such 
encounters are not very probable in stellar systems that evolve into 
star clusters like the Orion nebula cluster. This implies that star for- 
mation in UCDs may be influenced by processes that do not play 
a significant role in less massive stellar systems. The final masses 
and Plummer-radii resulting from the calculations in this paper are 
shown in Fig.l 161 

The possible initial conditions we uncover here (Figs. [74l 
and !15l > include densities as high as 10 8 Mq pc -3 for the forming 
UCDs with top-heavy IMFs (a ^ 1.9). The super nova rates are at 
times as high as one per year in the UCD-models with with an ini- 
tial stellar mass of 10 7 Mq and the most top-heavy IMFs (Fig. [4} 
and higher by a factor of 10 in the UCD-models with an initial stel- 
lar mass of 10 8 M Q . 

Starting from our initial conditions (Table [TJ, we seek those 
final models that represent UCDs in terms of their radii, masses 
and M/Ly -ratios in the following. 
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Figure 16. The final Plummer-radii against their final masses of all models 
that have not dissolved at the end of the integration. Different symbols en- 
code different IMFs, identified by the numbers assigned to them in Tab [2] 
The assumed star formation efficiency and heating efficiency is given at the 
top of each panel. The shaded regions indicate the parameter space occupied 
by real UCDs. 



• If the UCDs form as star clusters with a high star-formation 
efficiency and a high heating efficiency (as discussed with the case 
SFE=1 and HE=1; see Section [3~Tt . the properties of present-day 
UCDs are reproduced from models with all IMFs in Table ^ex- 
cept IMF 1, i.e. with stellar populations with high-mass IMFs in 
the whole range from a — 2.3 (canonical IMF) to a = 1.1 (see 
Tableland Fig. I16t. The different models imply however differ- 
ent ages and different stellar remnant populations for the UCDs, 
because they are constrained by the average M/Ly ratio that is 
observed for UCDs (cf. table 3 in DKB). A consistency check be- 
tween the models in this paper and the models in DKB is provided 
in Table|7]Note that the model from DKB where stars with an ini- 
tial mass larger than 25 Mq are assumed to evolve into black holes 
that have 50 per cent of the mass of their progenitors and all com- 
pact remnants are thought to be retained by the UCD is not listed 
in Table [7] This is because it is not consistent with the assumption 
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Table 7. Consistency-check between the models for the remnant populations of the UCDs discussed in DKB based on their dynamical M/Ly ratios with the 
models discussed here for their early dynamical evolution, assuming SFE=1 and HE=1 (for other SFEs and HEs see Section[4}- DKB consider two different 
ages for the UCDs. The first column specifies various remnant populations of the UCDs, as found in DKB. They differ by the mass of the SN remnants and 
which fraction of them remains bound to its host UCD. Concerning the kind of compact remnant a SN leaves, it is assumed that stars with initial masses 
between 8 Mq and 25 Mq become neutron stars with a mass of 1.35 Mq. Stars with initial masses above 25 Mq become black holes with either 10 or 50 per 
cent of the initial mass of their progenitors (t«bh = O.lor 0.5m*). The upper mass limit of the IMF is 100 Mg in all models. The second column displays 
the high-mass IMF slopes, 75, which correspond to these remnant populations, given the mean dynamical M/Ly ratio of the UCDs (see DKB for details). The 
uncertainties on a are calculated from the uncertainties on the mean M / Ly ratios. Columns 3 to 5 indicate how consistent the models discussed in DKB are 
with the ones discussed here. In this context, a '+' means that a ± 0.1 agrees with the high-mass IMF slope in one of the UCD-models marked with a '(+)' 
in Table|4] i.e. \a — a\ < 0.1, where a is the high-mass slope of the IMF whose number according to Table|2]is given in brackets. A 'O' has an analogous 
meaning, but it is only required that \a — a\ < 0.2. A ' — ' indicates that a model with these initial condidtions can only reproduce final parameters as observed 
in UCDs with an IMF with \a - a\ > 0.2. 
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that UCDs loose most of the mass that was initially locked up in 
their massive stars. However, it seems likely that UCDs loose in- 
deed most of this mass (see Section 12.2.21 The model from DKB 
where all matter is kept within the UCDs is also omitted from Ta- 
ble UJ even though it would seem consistent with the UCD-models 
with initial stellar mass of 10 s Mq, a SFE of 0.4 and a HE of 0.03. 
These initial parameters lead however to a unrealistic situation at 
the end of the calculation, because the UCD-models that can evolve 
this way stay too compact for being consistent with real UCDs (see 
SectionEHli. 

As a result of the comparison shown in TableUJ the models with 
IMFs 5 and 6 (canonical high-mass IMF slope) can be excluded as 
formation scenarios for the UCDs as a class of objects. This is be- 
cause these models suggest that the M/Ly ratio is consistent with 
the ones predicted by simple stellar population models, which is no t 
the case for UCDs toabringhausen et al.l2008l ; lMieske et alj2008h . 
The models with IMFs 5 and 6 would however be consistent with 
the UCDs in the Fornax Cluster if they are very old, because their 
average M/Ly ratio is som ewhat lower than the ones of UCDs in 
general jMieske et ai1l2008h . 

• If the UCDs form as star clusters with a moderate star- 
formation efficiency and high heating efficiency (as discussed with 
the case SFE=0.4 and HE=1 ; see Section [3T2l . extremely top-heavy 
IMFs (a = 1.1) can be excluded because they would lead to 
the complete dissolution of the cluster. Model m8_r5_hl_s04 with 
IMF 4 (a = 1.9) resembles a massive present-day UCD at the 
end of the integration. A comparison with table 3 in DKB shows 
that this model is consistent with two cases listed there. The first of 
them is the case of the UCDs being 13 Gyr old, keeping 20 per cent 
of the SN remnants and black holes, which retain 50 per cent of the 
mass of their progenitor stars. The second is the case of the UCDs 
being 13 Gyr old, keeping all SN remnants and black holes having 



10 per cent of the mass of their progenitor stars. These would be 
the only cases where a table analogous to Table [7] would indicate 
consistency between the UCD-models here and the ones in DKB. 
As in the case of SFE=1 and HE=1, the models with IMFs 5 and 6 
(canonical high-mass IMF slope) can be excluded as formation sce- 
narios for the UCDs as a class of objects, because of their to o-low 
M/Ly ratio jDabringhausen et all2008l ; lMieske et alj|2 008). 

• If the UCDs form as star clusters with a moderate star- 
formation efficiency and low heating efficiency (as discussed with 
the case SFE=0.4 and HE=0.03; see Section [3~3l >. the models with 
an initial mass of 2.5 x 10 s M lead to the unrealistic case that 
gas of the order of 10 8 Mq is confined on a very small volume at 
the end of our calculations. Models starting with an initial mass 
of 2.5 x 10 7 Mq on the other hand dissolve for the two most 
top-heavy IMFs, like in the case of SFE=0.4 and HE=1. Model 
m7 _r3_s04_h003 with IMF 4 (a = 1 .9) is similar to a small present- 
day UCD at the end of the calculation. A comparison with table 3 
in DKB shows that this model is consistent with two cases listed 
there. The first of them is the case of the UCDs being 13 Gyr old, 
keeping 20 per cent of the SN remnants and black holes, which re- 
tain 50 per cent of the mass of their progenitor stars. The second is 
the case of the UCDs being 13 Gyr old, keeping all SN remnants 
and black holes having 10 per cent of the mass of their progenitor 
stars. These would be the only cases where a table analogous to 
TableUJwould indicate consistency between the UCD-models here 
and the ones in DKB. 

Note the difference to the UCD-models with moderate SFE and 
high HE. In the case of a moderate SFE and high HE, consis- 
tency between the UCD-models from this paper and the mod- 
els for the remnant populations of UCDs from DKB is reached 
for a UCD-model starting with an initial stellar mass of 10 8 Mq, 
whereas in the case of a moderate SFE and a low HE consistency 
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is reached for a UCD-model starting with an initial stellar mass 
of 10 7 Mq. The UCD-models from this paper are however in both 
cases consistent with the same models for the remnant populations 
of UCDs in DKB. As in the case SFE=1 and HE=1, the models 
with IMFs 5 and 6 (canonical high-mass IMF slope) can be ex- 
cluded as formation scenarios for the UCDs as a class of objects , 
because of their to o-low M/Ly ratio iDab ringhause n"et alj|2008l ; 
iMieske et ai1l2008l) . 

Thus, in summary, the preferred solution of the initial condi- 
tions problem for a SFE of 0.4 are a proto-UCD with a stellar initial 
mass of 10 8 Mq and a projected half-mass radius of 5 pc (HE=1) or 
a proto-UCD with a stellar initial mass of 10 7 Mq and a projected 
half-mass radius of 3 pc (HE=0.03). UCD-models with SFE=1 are 
discussed in Section [3~T1 and for SFE=1 and HE=1, the preferred 
solutions are presented in Table [7] 

The comparison between the final parameters of our models 
and observed parameters of present-day UCDs contain some un- 
certainties for initial parameters that lead to the formation of UCDs 
because of a number of approximations and simplifying assump- 
tions (also see Section[2]for this matter). For instance, the d ensity 
profil es of UCDs are usually better fitted by a King profile J Kind 
1 19661) than by a Plummer profile. Also, the calculations performed 
here stop at 250 Myr, whereas UCDs are ~ 10 Gyr old. Thus UCDs 
will have suffered from adiabatic mass loss through the evolution 
of intermediate-mass stars, if the material expelled by them is not 
used up in the formation of subsequent stellar populations. Finally, 
the tidal field of the host galaxy of the UCD may play a role for its 
evolution on a Gyr time-scale. The performed comparison demon- 
strates however that also the more rapid early mass loss triggered 
by an over-abundance of massive stars does not necessarily lead to 
complete dissolution of massive, dense stellar systems, but can re- 
sult in objects similar to a UCD. The existence of UCDs is therefore 
not in contradiction with their formation with a top-heavy IMF. In 
a number of cases, a top-heavy IMF leads to a strong inflation of 
the modelled UCDs, but does not completely disintegrate them. 
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